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Water freezes into ice, atomic spins spontaneously align in a magnet, liquid helium becomes 
superfluid: Phase transitions are dramatic phenomena. However, despite the drastic change 
in the system's behaviour, observing the transition can sometimes be subtle. The hallmark 
of Bose-Einstein condensation (BEC) and superfluidity in trapped, weakly interacting Bose 
gases is the sudden appearance of a dense central core inside a thermal cloud 1-7 . In strongly 
interacting gases, such as the recently observed fermionic superfluids 8 , this clear separation 
between the superfluid and the normal parts of the cloud is no longer given. Condensates 
of fermion pairs could be detected only using magnetic field sweeps 9-11 into the weakly in- 
teracting regime. The quantitative description of these sweeps presents a major theoretical 
challenge. Here we demonstrate that the superfluid phase transition can be directly observed 
by sudden changes in the shape of the clouds, in complete analogy to the case of weakly in- 
teracting Bose gases. By preparing unequal mixtures of the two spin components involved in 
the pairing 12,13 , we greatly enhance the contrast between the superfluid core and the normal 
component. Furthermore, the non-interacting wings of excess atoms serve as a direct and 
reliable thermometer. Even in the normal state, strong interactions significantly deform the 
density profile of the majority spin component. We show that it is these interactions which 
drive the normal-to-superfluid transition at the critical population imbalance of 70(5)% 12 . 

The dramatic signature of BEC in weakly interacting gases in atom traps derives from a 
natural hierarchy of energy scales: The critical temperature for condensation T c oc n 2 / 3 at particle 
density n is much larger than the chemical potential (divided by the Boltzmann constant k B ) of 
a pure condensate, /J oc na, which measures the interaction strength between atoms (a is the 
scattering length). Hence, for weak (repulsive) interactions (na 3 <C 1), the condensate is clearly 
distinguished from the cloud of uncondensed atoms through its smaller size and higher density. 
However, as the interactions are increased, for example by tuning a using a Feshbach resonance, 
this hierarchy of energy scales breaks down, as fi can now become comparable to k B T c . In Fermi 
gases with weak attractive interaction (a < 0), /x oc Ep will even exceed ksTc oc Ep exp _7r//2fcF ' a ' 
by far (Ep = h 2 kp/2m is the Fermi energy). Both the normal and the condensed cloud will here 
be of the same size given by the Fermi Radius Rp oc \JE F . 

The phase transition from the normal to the superfluid state, although dramatic in its conse- 
quences, is thus not revealed by a major change in the appearance of the gas. Indeed, in strongly 
interacting Fermi gases no deviation from a normal cloud's shape could so far be detected; neither 
in the unitary regime, where a diverges, nor on the attractive (BCS-)side of a Feshbach resonance. 
Theoretical works predicted small "kinks" 14-16 or other slight deviations 17 in the density profiles 
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of the gas in the superfluid regime, but after line-of-sight integration these effects have so far been 
too small to be observable. Condensates could only be observed via rapid magnetic field ramps 
to the BEC-side (a > 0) of the Feshbach resonance, performed during expansion 91 °. This sud- 
denly reduced the condensate's chemical potential and let the thermal fraction grow beyond the 
condensate size. A similar ramp was used to detect vortices on resonance and on the BCS-side in 
the demonstration of fermionic superfluidity 8 . However, these magnetic field ramps are difficult 
to model theoretically, and a satisfactory quantitative comparison of e.g. the condensate fraction 
with experiments has not been accomplished 18 ~ 21 . 

In this work we demonstrate that the normal-to-superfiuid phase transition in a strongly in- 
teracting Fermi gas can be directly observed from absorption profiles, without the need of any 
magnetic field ramps. As in the case of weakly interacting BECs, preparation, expansion and de- 
tection of the sample all take place at the same, fixed magnetic field and scattering length. As for 
BECs, the phase transition is observed by a sudden change in the shape of the cloud during time-of- 
flight expansion, when the trap depth is decreased below a critical value. To clearly distinguish the 
superfluid from the normal component we break the number symmetry between spin up (majority 
atom number iVf) and spin down (minority atom number iVj_) and produce an unequal mixture of 
fermions (imbalance parameter 5 = (iVj — jVjJ/(jVf + iVj). Standard BCS superfluidity requires 
equal densities of the two spin components. Hence, when cooled below the phase transition the 
cloud should show a sudden onset of a superfluid region of equal densities. Indeed, below a critical 
temperature, we observe a bimodal density distribution of the smaller cloud. 

Breaking the symmetry in atom numbers thus produces a direct and striking signature of the 
superfluid phase transition. A similar situation has been encountered in Bose-Einstein condensa- 
tion, where breaking the symmetry of a spherical trap resulted in dramatic anisotropic expansion 
of the condensate, now a hallmark of the BEC phase transition. 

Fig. 1 shows column density profiles of the two imbalanced spin states for different points 
along the evaporation path corresponding to different temperatures, and for three magnetic fields 
corresponding to the BEC-side, exact resonance and the BCS-side of the resonance. For high final 
trap depths (upper panels in Fig. 1), the smaller cloud has the expected shape of a normal, non- 
superfluid gas: It is very well fit using a single, finite temperature Thomas-Fermi-profile (with 
central optical density, radius and the fugacity as independent fit-parameters). However, below a 
critical trap depth, a second, denser feature appears in the center of the minority component (lower 
panels in Fig. 1). This onset of bimodality occurs very suddenly as the trap depth is lowered, as can 
be seen from Fig. 2: Around the critical point, the atom number (Fig. 2a) and population imbalance 
(Fig. 2b) are practically constant, and the temperature (Fig. 2c) varies in a smooth linear way with 
the trap depth. In contrast, below the critical trap depth, the shape of the smaller cloud starts to 
deviate drastically from the Thomas-Fermi distribution of a normal gas, as quantified in Fig. 2d. 
This sudden increase in the standard deviation of a fit to a single-component fitting function is a 
standard way of identifying the BEC phase transition in a model-independent way 2 . 

Fig. 2e displays the fact that below the critical trap depth a new, third radius is required to 
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Figure 1: (Color online) Direct observation of the phase transition in a strongly interacting two- 
state mixture of fermions with imbalanced spin populations. Top a-c and bottom d-f rows show 
the normal and the superfluid state, respectively. Panels a and d were obtained in the BEC-regime 
(at 781 G), b,e on resonance (B = 834 G) and c,f on the BCS-side of the Feshbach resonance (at 
853 G). The profiles represent the azimuthal average of the column density after 10 ms (BEC- 
side) or 11 ms (resonance and BCS-side) of expansion. The appearance of a dense central fea- 
ture in the smaller component marks the onset of condensation. The condensate causes a clear 
depletion in the difference profiles (bottom of each panel). Both in the normal and in the super- 
fluid state, interactions between the two spin states are manifest in the strong deformation of the 
larger component. The dashed lines show Thomas-Fermi fits to the wings of the column den- 
sity. The radii and mark the Fermi radius of a ballistically expanding, non-interacting 
cloud with atom number iVj, N±. The trap depth U, the atom numbers, the population imbal- 
ance 5 = (iVj — A^)/(iV| + N±), the interaction parameter l/k F a, the temperature T and the 
reduced temperature T/T F were: a, U = 4.8 fjK, iV T = 1.8 x 10 7 , N t = 2.6 x 10 6 , 5 = 75%, 
l/k F a = 0.42, T = 350 nK, T/T F = 0.20. b, U = 3.2 ^K, iV T = 1.8 x 10 7 , N l = 4.2 x 10 6 , 
5 = 63%, l/k F a = (resonance), T = 260 nK, T/T F = 0.15. c,U = 2.5 fiK, iV T = 1.5 x 10 7 , 
N t = 4.5 x 10 6 , 5 = 52%, l/k F a = -0.13, T = 190 nK, T/T F = 0.12. d, U = 0.8 //K, 
iVf = 6.5 x 10 6 , N t = 1.5 x 10 6 , 5 = 62%, l/k F a = 0.67, T = 50 nK, T/T F < 0.05. e, 
U = 1.1 /iK, iV T = 1.5 x 10 7 , N l = 3.8 x 10 6 , S = 60%, l/k F a = (resonance), T = 70 nK, 
T/T F = 0.06. f , U = 1.2 fiK, = 1.3 x 10 7 , N t = AA x 10 6 , 5 = 50%, l/k F a = -0.15, 
T = 100 nK, T/T F = 0.08. 
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Figure 2: (Color online) Characterization of the phase transition. The data characterize the evolu- 
tion of the fermion mixture as the cloud is evaporatively cooled by lowering the trap depth. The 
chosen magnetic fields are identical to those in Fig. 1 . Data obtained from the majority (minority) 
cloud are shown as diamonds (circles). Shown are a) the atom number, b) population imbalance 
between the two spin states and c) the temperature of the spin mixture as determined from the 
non-interacting wings of the larger cloud's profile, d) A finite temperature Fermi-Dirac (for res- 
onance and the BCS-side) or gaussian (for the BEC-side) distribution is fit to the minority cloud. 
The phase transition is marked by a sudden increase in x 2 a s the condensate starts to appear, e) 
Outer radii of the majority and minority cloud (for the minority cloud on the BEC-side: thermal 
cloud radius, all other cases: Thomas-Fermi radius) as well as the condensate radius, defined as 
the position of the "kink" in the minority profile (see Fig. 1). The radii of the majority and mi- 
nority clouds are normalized to the Fermi-Radii ify j of non-interacting atoms with atom number 
iVj and adjusted for ballistic (hydrodynamic) expansion. Note that the imbalance decreases dur- 
ing evaporation because the larger majority cloud incurs stronger evaporative losses. For the data, 
three (BEC&Resonance) to five (BCS) independent measurements were averaged. 
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describe the two clouds. As we will see below, the appearance of this central feature coincides 
with the appearance of the fermion pair condensate in experiments involving the magnetic field 
ramp technique 1(M2 . It is this condensate which contains the superfluid vortices in 812 . We are 
thus naturally led to interpret the central core as the condensate of fermion pairs, and the outer 
wings as the normal, uncondensed part of the cloud. This constitutes the first direct observation of 
the normal-to-superfiuid phase transition in resonantly interacting Fermi gases on resonance and 
on the BCS-side (i.e. without a magnetic field sweep that so far cannot be quantitatively accounted 
for). 

Already at high temperatures, above the phase transition, the larger cloud's profile is strongly 
deformed in the presence of the smaller cloud, a direct signature of interaction. Indeed, on 
resonance the cloud size of the minority component is significantly smaller than that of a non- 
interacting sample with the same number of atoms (see Fig. 2e). At the phase transition, the outer 
radii of the clouds do not change abruptly. This demonstrates that interactions, not superfluidity, 
are the main mechanism behind the reduced cloud size of an interacting Fermi gas. 

On the BEC-side, the condensate is clearly visible in the larger cloud. On resonance, how- 
ever, the condensate is not easily discernible in the larger component's profiles at the scale of 
Fig. 1. Nevertheless, we have found a very faint but reproducible trace of the condensate when 
analyzing the curvature of these column density profiles (see Fig. SI in the Supplementary Infor- 
mation). On resonance and on the BCS-side, the onset of bimodality in the smaller cloud can be 
clearly observed for imbalances larger than ~ 20% (but below a certain critical imbalance, see 
below), for which the condensate is small compared to the minority cloud size. With increasing 
magnetic field on the BCS-side (i.e. with decreasing interaction strength), the bimodality becomes 
less pronounced and is not clearly discerned beyond 853 G (l/k F a < —0.15). 

Thermometry of strongly interacting Fermi gases has always been a major difficulty in ex- 
periments on strongly interacting fermions 22 . A thermometer can only be reliable if the working 
substance is not affected by the sample to be measured. In equal mixtures of fermions, the two 
overlapping atomic clouds are strongly interacting throughout. Temperatures determined from a 
non-interacting Thomas-Fermi fit to these clouds need calibration based on approximate theoretical 
calculations 22 . In addition, as will be reported elsewhere, we find that those fits do not describe the 
profiles of a partially superfluid Fermi gas as well as they do in the normal state, in agreement with 
theory 14-17 . In the case of imbalanced mixtures, the wings of the larger component, where the spin 
down species are absent, are non-interacting and thus serve as a direct thermometer (see Fig. 2c). 
For an imbalance of 5 = 75(3)%, we determine the critical temperature for the phase transition on 
the BEC-side at I /k F a = 0.46 to be T/T F = 0.18(3) (k B T F = ^(3(7V T + A^)) 1 / 3 - Fermi energy 
of a non-interacting, equal mixture with the same total number of fermions iVj + N±, oo - geometric 
mean of the trapping frequencies). This corresponds to T/T c ^ = 0.55(9) when comparing the 
temperature to the critical temperature T c> i for Bose-condensation in a non-interacting gas with 
Ni bosons. The reduction in the critical temperature is a direct consequence of strong repulsive 
interactions between the molecules. On resonance, at 5 = 59(3)%, we find T/T F = 0.12(2), and 
on the BCS-side l/k F a = -0.14 for 5 = 53(3)% we obtain T/T F = 0.11(2), where we have nor- 
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malized the temperature by the Fermi temperature of an equal mixture with the same number of 
atoms, k B T F = ftw(3(N^ + iVj) 1 / 3 {uj/2-k - geometric mean of the trapping frequencies). These 
are the first directly measured and reliable temperatures for the superfluid transition in strongly 
interacting Fermi gases. They may serve as a checkpoint for theoretical models. 

We note that the critical temperature will in general depend on the population imbalance. For 
example, for large enough imbalance on resonance or on the BCS-side, no condensate will form 
even at zero temperature 12 , as we discuss below. Here, the critical temperature for superfluidity 
will be zero. 

An important qualitative difference distinguishes the BEC-side from resonance: At the low- 
est temperatures on the BEC-side, the gas consists of only two parts: The superfluid core sur- 
rounded by a fully polarized degenerate Fermi gas of the excess species. On resonance and on the 
BCS-side, however, there exists a third region, a normal state in which both species are mixed. 
Several recent theories describe density profiles of imbalanced Fermi mixtures 23 ~ 28 . Mean-field 
theories that neglect interactions in the normal cloud and between the normal and condensed cloud, 
are only in qualitative agreement with our results. Descriptions which exclude the mixed region or 
find superfluidity on resonance at all population imbalances are ruled out by our observations. 

To elucidate the origin of the clear separation between condensate and normal components, 
we varied the population imbalance at our coldest temperatures and on resonance. Fig. 3b shows 
several resulting profiles after 1 1 ms expansion from the trap. For large imbalances, 5 > 70%, the 
minority cloud is not bimodal and well fit by a (unconstrained) Thomas-Fermi profile. At a critical 
imbalance of 5 70%, the condensate appears and then grows further as the imbalance is reduced 
(for the cloud radii see Fig. S2 in the Supplementary Information). 

To characterize the appearance of the condensate for imbalances around 5 = 70%, a Thomas- 
Fermi profile is fit to the wings of the minority cloud. The fraction of atoms not contained in this 
fit is a measure of the condensate fraction (see Fig. 3). We find a critical imbalance of S c = 
70(5)% above which the condensate disappears. This agrees with our previous work 12 , where we 
employed a rapid ramp method to the BEC-side to extract the condensate fraction. We observed the 
quantum phase transition from the superfluid to the normal state as a critical population imbalance 
of S c = 70% was exceeded. This strongly suggests that the bimodality observed here directly in 
the minority component and the bimodality observed in molecular clouds after a magnetic field 
sweep are signatures of the same phase transition. 

The transition at 5 C is known as the Clogston limit of superfluidity 12 ' 29 and occurs when the 
chemical potential difference S/i becomes larger than the (local) superfluid gap A(r) (see Supple- 
mentary Information). Here we present a simple picture for the character of this phase transition 
in a harmonic trap. Thomas-Fermi-fits for the normal clouds beyond 5 C allow a simple estimate 
of the central 3D-density of the gas (with an estimated accuracy of 20% for the relative density 
difference), shown in the inset of Fig. 3. For large imbalances, we find that the 3D-densities differ 
significantly, as is expected for two weakly interacting Fermi clouds. As the imbalance is re- 
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Figure 3: (Color online) Quantum phase transition to superfluidity for decreasing population imbal- 
ance, a, The "condensate fraction" of excess minority atoms, not contained in the Thomas-Fermi- 
fit, versus population imbalance on resonance, b, Column density profiles, azimuthally averaged, 
for varying population imbalance. The condensate is clearly visible in the minority component as 
the dense central feature on top of the normal background (finite-temperature Thomas-Fermi fit, 
dotted lines). Below the critical imbalance 5 C = 70%, the condensate starts to form. The inset in 
a shows the central densities of the larger (black diamonds) and smaller (grey circles) cloud in the 
normal state above 5 C . This demonstrates that here the central densities are unequal, suppressing 
superfluidity. The densities were calculated from the central optical density and the fitted size of 
the clouds, assuming local density approximation and adjusting for ballistic (hydrodynamic) ex- 
pansion of the outer radii of majority (minority) clouds. The data-points for the condensate fraction 
show the average of several independent measurements. 
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duced towards the critical 5 C , the central densities approach each other and become approximately 
equal around S c . This is a direct consequence of strong interactions in the normal state. In a non- 
interacting Fermi mixture with an imbalance of 5 C , the central densities would differ by a factor of 
2.4. 

This observation now offers an intriguing insight into the nature of a fermionic superfluid on 
resonance or on the BCS-side. Already in the normal state above Tc or beyond 5 = 5 C , interac- 
tions between the two spin states are strong. Indeed, this is directly seen in the deformation of the 
majority cloud due to the presence of the minority species (see Figs. 1, 3). However, here these 
interactions are not strong enough to let the central densities of the two clouds become compara- 
ble. At the critical imbalance the Clogston criterion Sfi = A(0) is fulfilled in the center of the 
trap. For smaller imbalance, a central superfluid region can form, the condensate. Its borders are 
defined by Sfi < A(r). The simple density estimate in Fig. |3] suggests that this region will be of 
equal densities, although more refined techniques to measure small density differences have to be 
developed to finally conclude on this question. Outside the superfluid region there is still a normal 
state with unequal densities of minority and majority components. The discontinuity in the clouds' 
densities at the normal-to- superfluid phase boundary gives rise to the visible kink in the column 
density profiles. Such a density discontinuity is characteristic for a first-order phase transition. 

Interestingly, most of the "work" to build the superfluid state is already done in the normal 
component by decreasing the density difference. Consequently, the critical population difference 
to form the superfluid is largely determined by the interactions in the normal gas. 

In conclusion, we have observed the normal-to-superfluid phase transition through the direct 
observation of condensation in an imbalanced Fermi mixture, on the BEC-side, on the BCS-side 
and right on the Feshbach resonance. Unequal mixtures offer a direct method of thermometry by 
analyzing the non-interacting wings of the majority species. Strong interactions are already vis- 
ible in the normal cloud as marked deformations of the majority profile. It is these interactions 
in the normal gas which squeeze the two components and eventually, at the critical imbalance, 
let them reach almost equal densities in the center, aiding the formation of the superfluid. Our 
method of direct detection of the condensate is a powerful new tool to characterize the superfluid 
phase transition. At the current level of precision, the appearance of a condensate after magnetic 
field sweeps and the direct observation of the central dense core occur together and indicate the 
normal-to-superfluid phase transition. An intriguing question is whether further phases are possi- 
ble, including a more exotic superfluid state with unequal densities. Several theories predict that 
the FFLO-state, a superfluid state with oscillating order parameter, should be present for imbal- 
anced spin populations 24 ' 30 . 

Methods 

Our experimental setup is described in previous publications 812 . A spin-polarized cloud of 6 Li 
fermions is cooled to degeneracy using a combination of laser cooling and sympathetic cooling 
with sodium atoms in a magnetic trap. After transfer into an optical trap, a variable spin mixture of 
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the lowest two hyperfine states, labelled |f) and ||) is prepared at a magnetic bias field of 875 G. 
Interactions between the two spin states can be freely tuned via a 300 G wide Feshbach resonance 
located at B = 834 G. At fields below B , two-body physics supports a stable molecular bound 
state (BEC-side), while at higher fields (BCS-side), no such bound state exists for two isolated 
atoms. Our trap combines a magnetic saddle potential with a weakly focused (waist w 120 jum) 
infrared laser beam (wavelength A = 1064 nm), leading to a harmonic axial confinement with 
oscillation frequency of v z = 22.8(0.2) Hz and a gaussian radial potential with variable trapping 
frequency v r in the central harmonic region. The trap depth U is related to v r and v z via 



The initial degeneracy of the spin-mixture is about T/T F ps 0.3. The strongly interacting gas 
is further cooled by decreasing the laser power of the optical trap in several seconds and evaporat- 
ing the most energetic particles. During the first few seconds, the magnetic field is adiabatically 
ramped to a chosen final field in the resonance region where the last stage of the evaporation 
(shown in Fig. 2) takes place. For detection, the optical trap is switched off and the gas expands in 
the remaining magnetic saddle point potential. After a variable time- of- flight an absorption image 
of atoms either in state |f) or |J,) is taken along the axial direction of the trap (the direction of 
the optical trapping beam). The cloud's radial symmetry allows for azimuthal averaging of the 
resulting column densities, leading to low-noise profiles 12 . 

For preparing clouds at the coldest temperatures (as shown in Fig. 3) with varying population 
imbalance, the spin mixture is evaporated down to a trap depth of 1 /iK over several seconds on 
resonance, after which the trap depth is increased again to 1.4 /xK for more harmonic confinement 
(trap frequencies: v r = 115(10) Hz and v z = 22.8(0.2) Hz). The temperature of the gas is 
determined to be T/T F < 0.06 for all 5 > 15%, and appeared to smoothly raise to T/T F = 0.11 
for an equal mixture, although thermometry in the interacting wings is problematic. The total atom 
number was 1.5 x 10 7 and constant to within 15% for all values of 5. 

The error in the critical temperature T C /T F for the phase transition is dominated by the 
uncertainty in the atom number entering the determination of T F , which we estimate to be 30% 12 . 
For Tp we use the harmonic approximation for the radially gaussian trapping potential, with the 
measured trapping frequencies reflecting the average curvature of the gaussian potential. The 
phase transition is observed above U = 2 /iK, where anharmonicities contribute only 3% to the 
error in T F . Note that anharmonicities do not affect the temperature measurement performed on the 
majority wings: Ballistic expansion of non-interacting atoms reveals their momentum distribution, 
regardless of the shape of the trap. 
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Supplementary Information 

Supplementary Figures 




Radius (mm) 



Supplementary Figure 4: (Color online) Signatures of the condensate on resonance in the spatial 
profiles. The curvature of the observed column density is encoded in shades of gray with white 
(black) corresponding to positive (negative) curvature. The outer radii of the two components and 
the condensate radius are shown as an overlay in the lower panel. As a direct consequence of 
strong interactions, the minority component causes a pronounced bulge in the majority density 
that is reflected in the rapid variation of the profile's curvature. The condensate is clearly visible in 
the minority component (5 > 0), but also leaves a faint trace in the majority component (5 < 0). 
The image was composed out of 216 individual azimuthally averaged column density profiles, 
smoothed to reduce technical noise. Data close to the cloud's center suffer from larger noise due 
to the lower number of averaged points. The central feature of about 50/im width is an artefact of 
smoothing in this region of increased noise. 
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Supplementary Figure 5: Outer radii of the two cloud profiles and condensate radius versus 
population imbalance. Data obtained from the majority (minority) cloud are shown as diamonds 
(circles). The outer radii of the clouds (black) are determined from Thomas-Fermi fits to the pro- 
files' wings, where the results of a zero-temperature and a finite temperature fit were averaged. For 
the minority cloud, the representative error bars indicate the difference between these two results. 
The position of the "bulge" in the majority profile (white diamonds) naturally follows the outer 
minority radius. The condensate radius is defined as the position of the "kink" in the minority 
profiles. It was obtained by a) fitting an increasing portion of the minority wings until a significant 
increase in \ 2 was observed (grey circles), and b) the position of the minimum in the profile's 
derivative (white circles). All sizes are scaled by the Fermi-radius of a non-interacting equal mix- 
ture. The minority radii were adjusted for the observed hydrodynamic expansion (expansion factor 
1 1.0). The non-interacting wings of the majority cloud expand ballistically (expansion factor 9.7), 
as long as they are found a factor 11/9.7 = 1.13 further out than the minority radius. For small 
imbalances (5 < 20%), also the majority wing's expansion will be affected by collisions. The grey 
diamonds give the majority cloud's outer radius if hydrodynamic expansion is assumed. 
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Supplementary Methods 
Hydrodynamic vs. ballistic expansion 

A non-interacting cloud of atoms simply expands ballistically from a trap. However, strongly 
interacting equal Fermi mixtures, above and below the phase transition, are collisionally dense and 
therefore expand according to hydrodynamic scaling laws 1-3 . These scaling laws only depend on 
the equation of state of the gas, e oc n 1 , with 7 = 1 for the BEC-side, 7 = 2/3 for resonance 
(a direct consequence of unitarity) and 7 = 2/3 for the BCS-side, away from resonance. In 
an unequal spin mixture of fermions, the expansion does not follow a simple scaling law. The 
minority cloud is always in contact with majority atoms and thus strongly interacting throughout 
the expansion, which is therefore hydrodynamic. The excess atoms in the wings of the larger 
cloud are non-interacting and will expand ballistically, as we have checked experimentally. The 
absorption images after expansion are taken along the axial direction of the trap (the direction of 
the optical trapping beam). In order to compare the expanded cloud sizes to the in-trap Fermi radii 
of non-interacting clouds (see Fig. 2 and Fig. S2 above) we scale the majority cloud with the 
ballistic factor for the radial direction 



where t is the expansion time and v z j \[2 gives the radial anti-trapping curvature of the magnetic 
saddle-point potential. The scaling factor for the hydrodynamic expansion of an equal mixture is 
given by the solution to a differential equation 2 3 . A priori, the minority cloud in unequal mixtures 
could expand with a different scaling, since the equation of state now depends on two densities. 
However, by imaging the cloud in trap and at different times during expansion, we found that 
the minority cloud's expansion is very well described by the scaling law for an equal mixture. 
In particular, the aspect ratio of the minority cloud did not change as a function of population 
imbalance (within our experimental error of 5%), and was equal to that of a balanced mixture. 

For the data on resonance in Figs. 3, SI and S2, which were obtained after 1 1 ms expansion 
out of a trap with radial (axial) frequency of v r = 113(10) Hz (y z = 22.8(0.2) Hz), the ballistic 
(hydrodynamic) expansion factor for the radial direction is 9.7 (1 1.0). 

Supplementary Discussion 
Signature of the condensate 

Fig. S 1 demonstrates that on resonance, the condensate is visible not only in the minority compo- 
nent, but also in the larger cloud as a small change in the profile's curvature. In the condensate 
region, the majority profile is slightly depleted when compared to the shape of a normal Fermi 
cloud. This effect is still significant on the BCS-side (see Fig. 1): Although here, the condensate 
is less visible in the smaller component than on resonance, the larger cloud's central depletion still 
produces a clear dip in the difference profile. 
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Radii in the unequal Fermi mixture 



Fig. S2 shows the outer radii of the majority and minority cloud, together with the condensate 
radius (on resonance, for the deepest evaporation compatible with constant total atom number ver- 
sus imbalance). As was the case for the phase transition at finite temperature, the outer cloud sizes 
change smoothly with imbalance. No drastic change is seen at the critical population imbalance. 
The radii are obtained by fitting the profiles' wings to the Thomas-Fermi expression for the radial 
column density n(r): 



with the central column density n , the fugacity A and the Thomas-Fermi radius R as the free 
parameters. Li 2 (x) is the Dilogarithm. The zero-temperature expression reduces to n(r) = n (l — 



Lower and upper bounds for the critical chemical potential difference at 5 C 

For the clouds at the critical imbalance 5 C , we now want to extract a lower and upper bound for 
the difference in chemical potentials 5/i c of the majority and minority component. This difference 
allows us to conclude that BCS-type superfluidity with imbalanced densities is not possible. 

The chemical potential difference 5fi = 2h = (/tf — /^) measures the energy cost, relative 
to pi = (/if + AH)/2, to add a particle to the cloud of excess fermions. A, the pairing gap, is the 
energy cost for this additional majority particle to enter the superfluid. Both the critical temperature 
T c and the critical chemical potential difference S/j, c provide a measure of the superfluid gap: 
The superfluid can be either destroyed by raising the temperature or by increasing the population 
imbalance. If h c = 5fi c /2 < A, excess atoms will always stay outside the superfluid, in the phase 
separated normal state. For h c > A, excess atoms can enter the superfluid for h c > h > A. 
Hence, superfluidity with unequal densities, if allowed via h c > A, would be favored at large 
population imbalance, contrary to the interpretation in 4 , where such a state was proposed for small 
population imbalance. A recent Monte-Carlo calculation 5 for the Clogston limit on resonance 
gives h c = 1.00(5)A = 0.50(5) E F and can thus not decide on the question of superfluidity with 
imbalanced densities. 

We can attempt to extract the chemical potential from the cloud sizes R^^ - taking into 
account hydrodynamic expansion for the minority cloud and ballistic expansion for the excess 
fermions. For the majority cloud, we find ji c ^ = l/2muj 2 R% = 1.21(6)E F . For the minority 
cloud, we find l/2mu; 2 Rl = 0.39(10)-Ep- Throughout the smaller cloud, minority atoms are 
always strongly attracted by majority atoms. This strong attractive interaction likely reduces their 
chemical potential from the above upper limit. The difference of the chemical potentials S/j, c = 2h c 
is thus given by h c = (fj, c ^ — /i Cj jJ/2 > 0A1(6)E F = 0.51//, our lower bound. Another condition 
on h c concerns whether the normal state can be mixed, h c < /t, (minority and majority atoms in 
the same spatial region) or whether the normal state is always completely polarized h c > /t. Our 




r 2 /R 2 ) 2 . 
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observation of the mixed region in Fig. 1 immediately results in h c < /x, the upper bound. 

On resonance, A = 1.16/x in BCS-theory, while a recent Monte-Carlo study 5 obtains 
A = 1.2/i. If A > fi holds true, our finding of the upper bound on h c would imply h c < A 
and hence would exclude a superfluid with unequal spin densities (at least on the basis of BCS- 
theory, see 6 for a recent suggestion which goes beyond BCS). 
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